On the Properties of Plastic Ablators in Laser-Driven Material Dynamics Experiments 
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Radiation liydrodynamics simulations were used to study the effect of plastic ablators in laser- 
driven shock experiments. The sensitivity to composition and equation of state was found to be 
5-10% in ablation pressure. As was found for metals, a laser pulse of constant irradiance gave a 
pressure history which decreased by several percent per nanosecond. The pressure history could be 
made more constant by adjusting the irradiance history. The impedance mismatch with the sample 
gave an increase o(100%) in the pressure transmitted into the sample, for a reduction of several 
tens of percent in the duration of the peak load applied to the sample, and structured the release 
history by adding a release step to a pressure close to the ablation pressure. Algebraic relations 
were found between the laser pulse duration, the ablator thickness, and the duration of the peak 
pressure applied to the sample, involving quantities calculated from the equations of state of the 
ablator and sample using shock dynamics. 
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I. INTRODUCTION 



Studies of the properties and response of matter under 
extreme conditions rely increasingly on lasers to induce 
dynamic loading by ablation. Compared with the canon- 
ical technique of impact-induced shocks, laser ablation is 
more convenient for studying phenomena on shorter time 
scales (nanoseconds rather than microseconds), for inves- 
tigating the properties of single crystals and the basic mi- 
crostructural processes, and for developing new diagnos- 
tics. However, laser ablation has potential disadvantages 
including uncertainties in loading history, spatial varia- 
tions in irradiance from speckles or the overall intensity 
envelope, and the possibility of heating the sample under 
investigation by hot electrons or Bremsstrahlung x-rays. 
A thorough understanding of the processes of ablative 
loading is necessary to assure that the response of the 
sample can be distinguished from phenomena associated 
with the ablation process. 

In many experiments, dynamic loading can be induced 
by ablation of the sample material itself. Previous stud- 
ies have demonstrated that spatial variations and pre- 
heat were negligible, for ablatively-generated shocks from 
^100 MPa to ITPa in elements [l| and compounds 2], 
and for ramp compression In other experiments, it is 
highly desirable to use the laser energy to ablate a well- 
controlled layer of a different material, in contact with 
the sample. Advantages of a separate ablator include 
the ability to restrict multi-electron-volt ablation tem- 
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peratures to elements of low atomic number irrespective 
of the sample material, reducing the risk of generating 
hard x-rays; smoothing of small-scale spatial variations 
by beam dispersion, radiation transport, and mechanical 
equilibration within the ablation plasma; and the abil- 
ity to use a single ablator of known plasma properties 
for experiments on less well-understood sample materi- 
als. A current interest in dynamic loading experiments is 
the use of x-ray diffraction to study the sample material 
in the compressed state With ramp compression, the 
material of highest compression is closest to the ablation 
surface, so it is desirable to pass the x-rays in and out 
through the ablator and the ablation plume. It is thus 
particularly desirable to restrict the ablator to elements 
of low atomic number, dissimilar to the sample. When 
the sample material is ablated directly, heat conduction 
from the plasma and re-condensation of ablated material 
at the end of the laser pulse can lead to material with 
a locally different microstructure @, an added compli- 
cation when interpreting the microstructure of recovered 
samples. A dissimilar ablator acts as a palliative in this 
case. By appropriate choice of material, an ablator may 
also increase the efficiency of conversion of laser energy 
to pressure in the sample. 

An ablator layer can however complicate dynamic load- 
ing experiments. The time at which the load is first 
applied to the sample depends on the ablator thickness 
and the laser irradiance. The mechanical impedance mis- 
match between the ablator and the sample causes wave 
interactions, leading to multiple waves in the sample 
which must be discriminated from the response of the 
sample. If the composition and properties of the ablator 
material are not adequately known, it introduces addi- 
tional uncertainty in the experiments. 
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Plastics are attractive ablator materials as they are 
stable to handle and readily available with compo- 
sitions of low atomic number, such as polystyrene 
(C8H8)„ and polyparaxylylene-N (parylene-N) (C8H8)„, 
and polyethene (C2H4)„. Deposition and coating tech- 
niques exist which can produce layers of uniform and 
well-characterized thickness 0|- Polystyrene and 
parylene-N ablators have been used as a coating for glass 
spheres in inertial confinement fusion experiments. Pre- 
vious studies have investigated the properties of the ab- 
lation plasma including two dimensional effects in pla- 
nar targets [1] , and have investigated the generation and 
mitigation of preheat in the sample ■ Significant 

discrepancies have been identified between pressures in- 
ferred from experiments and predicted in radiation hy- 
drodynamics simulations Parylene layers have been 
used in material dynamics experiments at lower irradi- 
ances, again with some discrepancy between simulation 
and experiment [T]| . 

Here we present studies of the behavior of plastic abla- 
tors under conditions relevant to shock physics and mate- 
rial dynamics experiments, including sensitivities to un- 
certainties in the composition or equation of state of the 
ablator, and wave interactions caused by the impedance 
mismatch with the sample. The effect of wave interac- 
tions is presented in a form convenient for the design of 
laser-driven shock experiments. 



II. RADIATION HYDRODYNAMICS 
SIMULATIONS 

The response of a plastic layer to intense laser irradi- 
ance, and the effect of the impedance mismatch between 
the plastic and a denser sample material behind, were in- 
vestigated using numerical simulations. Laser-matter in- 
teraction is non-linear and time-dependent, so spatially- 
resolved radiation hydrodynamics is needed to integrate 
the coupled radiation transport and continuum dynam- 
ics equations forward in time. In the regime of interest, 
laser ablation and dynamic loading can be simulated ac- 
curately by assuming three temperature hydrodynamics 
(ions, electrons, and radiation), including thermal con- 
duction and radiation diffusion, and calculating the ab- 
sorption of the laser energy in the expanding plasma 
cloud through the electrical conductivity Thermal 
conduction and radiation diffusion are necessary only 
near the ablated surface of the plastic: the bulk of the 
ablator and the sample behind are affected only by the 
hydrodynamics. 

Material dynamics experiments are designed to apply 
a simple loading history to the sample: ideally close to 
one-dimensional in space, i.e. with a time-dependence 
that is independent of position over a region of the sam- 
ple surface. Experiments are usually designed so that 
edge effects, two- or three-dimensional dependencies in 
the loading history, do not affect some portion of the 
sample for long enough for the one-dimensional response 



to be measured. We therefore use one-dimensional radi- 
ation hydrodynamics, and estimate the size of the region 
affected by the edges. 

Radiation hydrodynamics is the simultaneous solution 
of the continuum dynamics, heat conduction, and radia- 
tion transport equations [12, [l^l ■ Using the Lagrangian 
frame of reference, in which time derivatives are ex- 
pressed along characteristics moving locally with the ma- 
terial, the equations of continuum dynamics are 
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where t is time, p the mass density, u the particle velocity, 
p the pressure, and e the specific internal energy. These 
equations are closed by a mechanical equation of state 
(EOS) p{p,e). Each equation may also have a source 
term; here, heat conduction and radiation transport pro- 
vide source terms in the velocity and energy equations. 
Heat conduction and radiation transport involve temper- 
atures for the ions T^, electrons Tg, and radiation field T^. 
In the regime of interest, laser and thermal radiation cou- 
ples principally to the electrons, which equilibrate with 
the ions according to a compression- and temperature- 
dependent rate Tei{p,T). A thermal EOS relates e to 
Tf, in practice we used an EOS of the form {p, e}(p, T^), 
using Ti as one of the state parameters and calculating 
e(T) and Ti(e) when needed in the hydrodynamics. The 
equation of heat conduction is 



dT 

— = div [k{p, r,)gradr] , 



(4) 



where K{p,Ti) is the thermal conductivity. Radiation 
transport was calculated for an equilibrium radiation dis- 
tribution (represented by a single radiation temperature 
field Tr) The equation of single-temperature radiation 
diffusion is 



dU 

— = -prdivu 
ot 



divF -I- AnaaBp — c{aa 



)U, (5) 



where U = AaT^/c, eta and cts are frequency-averaged 
absorption and scattering coefficients. Bp is the inte- 
grated power in the Planck distribution at T^, the ra- 
diation flux 
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a is the Stefan-Boltzmann constant, and c is the speed of 
light. The absorption coefficient was approximated as the 
Rosseland mean opacity, and the emission coefficient as 
the Planck mean opacity. The transport and deposition 
of incident laser energy was calculated as the attenuation 
of a free-propagating flux, using the electrical conductiv- 
ity to calculate the attenuation, with energy deposited as 
heating of the electrons. 
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In our context, radiation hydrodynamics is an ini- 
tial value problem, where the material fields p and e 
are specified over a region R at some time to, time- 
dependent boundary conditions are specified for the con- 
tinuum PdR{t) or UQnit) and radiation flux Idnit), and 
the continuum and radiation fields are integrated for 
t > to. 

Simulations were performed using the HYADES radi- 
ation hydrocode, version 01.05.11 14]. This hydrocode 
used a one-dimensional (ID) Lagrangian finite-difference 
discretization of the material and leapfrog time integra- 
tion, with shock waves stabilized using artificial viscos- 
ity. The radiation fiux limiter was set to 0.03 of the free 
stream value - a common choice for simulations of this 
type [3. 



III. PROPERTIES OF PLASTIC ABLATORS 

The ablation properties of plastic coatings were inves- 
tigated by radiation hydrodynamics simulations as de- 
scribed above, to predict the relationship between irradi- 
ance and pressure, ablation rates, and the sensitivity to 
composition of the plastic. 

The properties of the plastic were represented through 
models for its pressure- volume-energy equation of state 
(EOS), its conductivities, and its opacity. Wide-ranging 
EOS were taken from the SESAME database [1^. 
Opacities for radiation diffusion were also taken from 
SESAME. Conductivities for laser deposition and heat 
conduction were calculated using the Thomas- Fermi ion- 
ization model ,14.,17j: this was found previously to be ac- 
curate for direct drive shock simulations on samples of a 
wide range of atomic numbers [l[ . A significant difficulty 
with parylene-N is the only EOS model available was con- 
structed by density scaling from an EOS for parylenc-C, 
which was fitted to shock Hugoniot data at 300-600 GPa 
with an uncertainty ~10% [l3,[l^. In parylene-C, one H 
is substituted by CI, and density scaling will not neces- 
sarily capture the different chemical behavior. To inves- 
tigate the uncertainty in pressures predicted in ablative 
loading experiments, the sensitivity to EOS and opacity 
was investigated by performing simulations for parylene- 
C and parylene-D, in which an additional H atom is 
substituted by CI, and also for polystyrene, which has 
the same composition as parylene-N though a different 
molecular structure. The molecular structure should be- 
come unimportant once the plastic has been converted 
to plasma. The models for parylene-D (SESAME tables 
7770 and 17770) [l^l were developed earlier and have 
been used more commonly in radiation hydrodynamics 
simulations, so these models were used for the studies of 
general ablation behavior described next. 

To allow a fine spatial resolution of the surface to be 
ablated, without requiring the same resolution in the bulk 
of the plastic where it remains solid, the material was 
discretized with geometrically-expanding cells. The ex- 
pansion factor between successive cells was around 1-5%: 



small enough to avoid inaccuracy in the hydrodynamic 
equations. The smallest cell, at the surface where abla- 
tion starts, was chosen to have an initial size in the few 
nanometer range. Numerical convergence with respect to 
spatial resolution was tested by performing simulations 
with 200 or 300 zones, with an expansion ratio of 2.5% 
or 2.0% respectively, giving zone widths of around 9.0 
or 2.6 nm at the ablation surface. The ablation pressure 
was around 5% higher in the finer-mesh simulations, indi- 
cating an adequate level of mesh convergence. The finer 
mesh was used for all subsequent simulations. 

For the ablation rate and pressure study, the plastic 
was chosen to be 50 /im thick, and the laser pulse was 
10 ns long with constant irradiance. To avoid numerical 
oscillations, the irradiance was taken to rise from and fall 
to zero over 0.1ns, which is typical of the laser systems 
used for the material dynamics experiments of interest 
here. The critical surface of the plasma was predicted to 
remain within a few micrometers of the initial position 
of the ablator, so the plasma plume should be essentially 
ID. Release and cooling at late times would be faster 
than predicted using ID simulations, because of lateral 
expansion. 

In the simulations, the surface of the plastic was heated 
to plasma conditions rapidly after the start of the laser 
pulse. The plasma then absorbed the incident laser en- 
ergy, forming a nearly steady state plume close to the 
plastic. The plasma pressure and reaction to plasma flow 
induced an elevated pressure in the remaining condensed- 
phase plastic, driving a shock wave. As was found previ- 
ously [H, [3| J a laser pulse of constant irradiance did not 
induce a constant driving pressure in the condensed re- 
gion. Instead, there was an initial spike of pressure tens 
of percent higher, decaying rapidly to a more constant 
sustained pressure. This sustained pressure typically de- 
creased by a few percent during the laser pulse. The EOS 
used exhibited simple, concave behavior, so the decaying 
pressure in the spike and the plateau was able to catch up 
and erode the higher pressure states ahead, reducing the 
shock pressure monotonically as it propagated through 
the material. [27j After the end of the pulse, the ablation 
pressure fell rapidly, and a rarefaction wave propagated 
into the shocked material, releasing the pressure. For 
the thickness and pulse duration chosen, and over the 
irradiance range considered, the shock reached the rear 
surface of the plastic before the rarefaction caught up, 
so a second rarefaction wave from the rear surface also 
propagated into the material. (Fig. [T|) 

Temperatures in the ablation plume were predicted 
to be of order 10-100 eV. The mass density and time 
scales were low enough that conductive heating of the 
condensed-phase ablator was not significant. (Fig. [2|) 

The relationship between irradiance and ablation pres- 
sure was investigated by performing simulations with the 
same 10 ns laser pulse duration, and varying the irra- 
diance between 0.1 and lOOPW/m^.f!!] The sustained 
pressure p following the initial peak was represented well 
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FIG. 1: Predicted pressure profiles in parylene-N at intervals 
of 1 ns, for an irradiance of 0.1 PW/m^ applied for 10 ns. The 
undisturbed material extends from to 50 /im, and the laser 
drive impinges from the left. 
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FIG. 2: Predicted profiles of ion temperature in parylene-N 
at intervals of Ins, for an irradiance of O.lPW/m^ applied 
for 10 ns. The undisturbed material extends from to 50 /im, 
and the laser drive impinges from the left. 



by a power law in irradiance /: 

p(GPa) = 8.614 [/(PW/m^)] 



0.833 



(7) 



The fitting uncertainty in the parameters is 1.5%, and the 
uncertainty in the irradiance-pressure values is around 
5%, because of temporal variations in the pressure. This 
result is consistent with previous studies on elemental 
metals Jj, though the exponent here is larger, and the 
irradiance-pressure relation is better reproduced by the 
power law fit despite covering a wider range of irradi- 
ance. The exponent is consistent with experimental and 
calculational results obtained for polystyrene Q , and the 



prefactor is similar to the value implied by these results. 
The exponent is also consistent with more recent simu- 
lation results for parylene [l3| . 

In this regime of irradiance, laser and thermal radia- 
tion penetrates the material much more slowly than the 
shock wave passes through, so there is a region of com- 
pressed material heated only by the passage of the shock 
between the ablation surface and the undisturbed mate- 
rial ahead of the shock. It is important to know the rate 
at which material is ablated, to ensure that the ablator 
is thick enough for a laser pulse of given irradiance and 
duration. The ablation rate was extracted from the sim- 
ulations by observing the time of first outward motion 
of points in the material. For convenience in choosing 
a thickness of ablator, a Lagrangian ablation rate was 
calculated, i.e. with respect to the original position of 
the material. The ablation-supported shock compresses 
material before it is ablated, so the ablation rate, with 
respect to the instantaneously compressed material, is 
lower than the Lagrangian rate. For laser irradiances 
from 0.1 and lOPW/m^, the ablation rates were fitted 
well by a straight line, indicating that a constant irradi- 
ance produces a constant ablation rate. At an irradiance 
of lOOPW/m^, the ablation rate decreased significantly 
over the 10 ns of the laser pulse. Ablation was predicted 
to continue at a similar rate for around 20% of the dura- 
tion of the laser pulse after its end, because of retained 
heat in the ablation plume. After this time, the ablation 
rate decreased rapidly (Fig. [3]) . The individual ablation 
rates ua were reproduced accurately by the relation 



UA(km/s) = 0.108 [/(PW/m^)] 



0.651 



(8) 



vahd for / between 0.1 and lOOPW/m^. Above 
lOPW/m^, this relation ceases to be valid over the full 
10 ns investigated: at lOOPW/m^ it is valid for the first 
Ins. The fitting uncertainty in the parameters is 0.5%. 
Note that 1 km/s=l /im/ns. This result is similar to 
the results reported for the ablation of polystyrene [1], 
though our simulations predict a larger prefactor. In 
the previous results, the depth of material ablated was 
inferred indirectly from measurements of the ablation 
plasma. [1^ 

There is significant uncertainty in the EOS of matter, 
particularly in the warm dense matter regime of the ab- 
lation plume. The uncertainty is correspondingly greater 
for a material for which no specific EOS has been devel- 
oped, as discussed above. It was found previously [ll, 0| 
that pressures induced in the condensed region were fairly 
insensitive to the details of the plasma EOS. The sensi- 
tivity was investigated here by performing equivalent sets 
of simulations using the EOS and opacity for parylene- 
C from the SESAME hbrary: (tables 7771 and 17620) 
[131, and also polystyrene (tables 7590 and 17593) [2l|. 
The two parylene models gave the same ablation pres- 
sure to within 1-2%, and the shock speeds were almost 
identical. The polystyrene model gave a smaller initial 
pressure spike and a sustained pressure that was up to 
several percent lower, but with less of a decrease over the 
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FIG. 3: Penetration of ablation surface during the laser pulse, 
for different irradiances. 
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FIG. 4: Sustained ablation pressure predicted for different 
plastics, with respect to the pressures calculated for parylene- 
D or polystyrene. 



pulse duration and hence a slower decay as it propagated 
through the plastic. The shock speed in polystyrene was 
significantly different: slower at low pressures and faster 
at high pressures. These variations are all small com- 
pared with typical uncertainties in mean irradiance in 
laser ablation experiments. (Fig. 21) 

The sensitivity to ablator material was investi- 
gated further by performing equivalent simulations for 
polyethene (tables 7171 and 17171) ^22]. The ablation 
pressures were 5-10% lower than for parylene-D, but 
within 5% of polystyrene. The pressure profile was gen- 
erally similar to that of polystyrene. (Fig.[4l) 

For x-ray diffraction experiments, the timing of the x- 
ray pulse with respect to the shock-loading of the sample 
is particularly important. The effect of composition and 
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FIG. 5: Shock speed in the ablator as a function of sustained 
ablation pressure, for different ablator compositions and equa- 
tions of state. For polyethene, the solid and dashed lines are 
Steinberg's analytical equation of state and SESAME table 
7171 respectively. For parylene, solid and dashed lines are for 
parylene-C and D respectively. 



EOS was significantly more sensitive than the ablation 
pressure (Fig. [5]). 

Experimentally, it has been found that some laser light 
may be transmitted through the plastic coating before 
the opaque plasma sheath forms, inducing earlier abla- 
tion of the sample material behind by a small amount 
of the laser energy. This early transmission was not pre- 
dicted in the radiation hydrodynamics simulations, which 
(as is common) did not include an adequate model of the 
line opacity and breakdown of the plastic. A palliative 
is to deposit a thin, metal flashing on the surface of the 
plastic to prevent early transmission. A typical flashing 
is 100 nm of Al, though other metals have been used 
The Al layer potentially has an effect on the ablation 
pressure, apart from the palliative effect. This effect was 
predicted by performing simulations in which the first 
100 nm of the plastic was replaced by Al. At O.lPW/m^, 
the presence of the Al increased the sustained ablation 
pressure by 10-15% and made it more constant, but in- 
troduced slight wave reverberations in the compressed 
ablator (Fig. [5]). The difference was smaller at 1 PW/m^, 
and negligible at lOPW/m^ and above. 

There is also a potential effect from x-ray preheating 
of the sample, induced by K-shell radiation from the ab- 
lator or flashing. X-ray preheating is very sensitive to 
intensity variations in the laser beam, as x-ray yield is 
highly non- linear in laser irradiance [23|. The fluence of 
preheating x-rays is therefore very sensitive to the beam 
quality of the specific laser system, including any phase 
plate used to reduce large-scale intensity variations. The 
degree of x-ray heating also depends on the opacity of the 
sample material. Experimental evidence and theoretical 
predictions for x-ray preheating in such experiments is 
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FIG. 6: Predicted pressure profiles in parylene-N at intervals 
of 2 ns, with (solid lines) and without (dashed lines) a flashing 
of 100 nm Al on the surface of the plastic. The laser irradiance 
was O.lPW/m^ applied for 10 ns. The undisturbed material 
extends from to 50 /im, and the laser drive impinges from 
the left. 

mixed: simulations have indicated that preheating; may 
be significant in our regime of interest [lO| , but whereas 
experiments have shown little or no preheating [H, 0, [E, [1] . 
While preheating is a potential concern when the initial 
temperature of the sample is important, it has a much 
smaller effect on wave interactions caused by the pres- 
ence of the ablator, so we do not consider it further in 
the present work. 

At the end of the laser pulse, the pressure at the abla- 
tion surface drops rapidly as the ablation plume expands 
and cools. Until two-dimensional expansion takes effect, 
the pressure gradient in the ablation plume is inversely 
proportional to the pulse length, so the rate at which 
the pressure falls is inversely proportional to the pulse 
length. 



FIG. 7: Example laser irradiance history designed to give a 
more constant shock pressure than from a constant irradiance 
of O.lPW/m^ 

pulse where the time-dependent irradiance was repre- 
sented in tabular form, with table entries at the start 
and end of the pulse of equal value for a constant irra- 
diance. The procedure followed was first to adjust the 
irradiance at the end of the pulse, giving a linear ramp 
starting at the initial irradiance, and then to add extra 
points between the first and last and change their val- 
ues to remove shorter time scale variations in pressure. 
The sustained ablation pressure was made constant by 
ramping the irradiance up in time, and the initial pres- 
sure spike was removed by reducing the irradiance early 
in time and following a more gentle approach to the main 
part of the pulse. This procedure was demonstrated pre- 
viously for direct ablation of Be ^23| . Calculationally, it 
was possible to make the induced pressure constant to 
any desired precision by adjusting the laser irradiance 
history (Fig. [7|). For irradiances from 0.1 to lOPW/m^, 
the ramp was 10-30% of the mean irradiance. 



IV. LASER PULSE SHAPING FOR CONSTANT 
ABLATION PRESSURE 

The ablation pressure can be made constant in time by 
adjusting the temporal shape of the laser pulse. Many, 
though not all, large laser systems allow the power his- 
tory to be controlled, so a desired temporal shape can be 
delivered to some finite level of precision. On lasers with 
large numbers of beams but a simple shape from each 
(such as a constant power), a smooth temporal shape 
can be approximated by altering the energy and relative 
timing of the beams. 

The principle of pulse shaping for constant ablation 
pressure was demonstrated for parylene-N ablators by 
performing simulations in which the laser irradiance was 
adjusted. These simulations were performed for a laser 



V. WAVE INTERACTIONS WITH THE 
SAMPLE 

When the ablator is in contact with a sample, the 
impedance mismatch between the ablator and the sam- 
ple lead to wave interactions which change the pressure 
history in the ablator and may lead to multiple waves in 
the sample. In most cases of interest, the sample has a 
higher shock impedance than the plastic ablator, so the 
ablation-induced shock reflects a shock from the interface 
between the ablator and the sample. 

In this discussion of wave interactions, we consider a 
simplified situation where the laser irradiance history has 
been adjusted to give a constant drive pressure in the ab- 
lator, which will induce a constant initial shock pressure 
in the sample. Early in the laser pulse, an ablation-driven 
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FIG. 8: Schematic of wave interactions in position-time space 
for laser loading in one dimension with an ablator of lower 
impedance than the sample. Diagram shows the ablation- 
driven shock (AS), transmitted shock (TS) in the sample, 
reflected shock (RS) in the ablator; ablation surface release 
(ASR), transmitted release (TR), reflected release (RR), and 
free surface release (FSR); ablation release (AR); ablation 
front (AF). Darker shades represent material of higher mass 
density. For simplicity, release waves are shown as lines rather 
than fans. 



shock is induced, driving a transmitted shock in the sam- 
ple and a reflected shock in the ablator. If the laser pulse 
is long enough, when the reflected shock reaches the ab- 
lation surface, double-shocked ablator material is then 
ablated, resulting in an ablation pressure that is slightly 
higher than before, but generally lower than the double- 
shocked pressure. The resulting ablation surface release 
wave propagates through the ablator and is transmitted 
into the sample, reflecting a further weak release wave 
back through the ablator. When the transmitted shock 
reaches the free surface of the sample, a release wave 
propagates backward into the sample. When the laser 
pulse ends, a strong release wave propagates through any 
remaining ablator and forward into the sample. Where 
the strong release waves interact, the sample material is 
subjected to tension, and spall may occur. (Figs [5] and 

E) 

The impedance mismatch between most metals and 
plastics such as parylene and polyethene means that a 
given ablation pressure in the plastic induces a consider- 
ably higher pressure in a metal sample. The pressure en- 
hancement was estimated for the plastic EOS used above 
and published EOS for Al and Cu, by finding the inter- 
section between the secondary shock Hugoniot for a given 
initial shock state in the plastic and the principal Hugo- 
niot in the metal in pressure-particle speed space. The 
EOS of Al was as above; that for Cu was SESAME ta- 



FIG. 9: Schematic of wave interactions in particle speed- 
normal stress space for laser loading in one dimension with 
an ablator of lower impedance than the sample. Loci for the 
ablator (A) and sample (S) are solid and dashed respectively. 
The principal shock Hugoniots are denoted HI. Ablation in- 
duces a state apl on the principal Hugoniot of the ablator. 
Because of the higher impedance of the sample, the state in- 
duced is the intersection of the sample's principal Hugoniot 
with the secondary Hugoniot H2 of the ablator. When the 
reflected shock reaches the ablation surface, a release wave 
R2 propagates through the ablator, releasing to the ablation 
pressure of the re-shocked ablator material, ap2. When this 
release wave reaches the sample, the ablator material must 
release further (R3) to intersect the release adiabat R2 of the 
sample. This intersection gives the next state induced in the 
sample. 



ble 3330. Shock Hugoniots and their intersections were 
found numerically [2a]. Repeating this procedure for a 
range of ablation shock pressures, relations were found 
between the pressure in the plastic and the pressure in 
the metal (Fig. fTO)) . 

Radiation hydrodynamics simulations were performed 
to predict the integrated effect of ablative loading via a 
plastic ablator on the loading history experienced by the 
sample. As for the ablator-only simulations, a geometri- 
cally expanding spatial discretization was used. The sam- 
ple was chosen to be 25 /im thick. Simulations were per- 
formed for ablators 10 and 20 /xm thick, with Al flashing. 
The laser pulse was taken to be 3 ns long, with constant 
irradiance, which are specifications often used for mate- 
rial dynamics experiments at the Janus laser [llj. With 
a laser pulse of constant irradiance, the pressure was pre- 
dicted to decrease with time during the pulse, though the 
compression varied less. Simulations were performed for 
irradiances of 5, 10, 20, and 50PW/m^. The peak pres- 
sure in the sample was in agreement with the impedance 
mismatch calculations, and decayed as the shock propa- 
gated through the sample. With a thicker ablator, the 
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FIG. 10: Effect of shock impedance mismatch for plastic ab- 
lators driving shocks into AI and Cu samples. The sample 
pressure ratio is the pressure in the sample divided by the 
ablation pressure. 



pressure in the sample was more constant. (Figs [TT] to 

For comparison, simulations of ablative loading of un- 
coated samples were performed. For a given irradiance, 
the pressure induced in the sample was significantly lower 
when ablated directly, though the duration of the period 
of high pressure was greater. Ablation of plastic was pre- 
dicted to generate a slightly lower pressure than ablation 
of the metal sample, but the large impedance mismatch 
more than compensated for the reduced ablation pres- 
sure. With an ablator, the pressure fell more slowly at 
the end of the laser pulse. With an ablator, the pres- 
sure varied less during the initial shock, with a constant 
laser irradiance. The inclusion of the low density abla- 
tor increased the coupling of laser energy to the sample. 
(Figs [Hand [HI) 

One way to think of the function of the ablator is that 
it converts laser energy into kinetic and potential en- 
ergy in the moving, compressed ablator, which is then 
partly transferred to the sample. For a given irradiance, 
laser pulse length, and sample material, there is an op- 
timum thickness for any ablator to maximize the dura- 
tion of the peak shock pressure applied to the sample. 
The optimum thickness is such that the laser pulse ends 
when the reflected shock reaches the ablation surface. 
When the ablation rate is a constant in time, this rela- 
tionship can usefully be expressed as an effective speed: 
microns of ablator per nanoseconds of laser pulse; this is 
a useful quantity when designing ablation experiments. 
The optimum relationship between thickness and pulse 
length depends on the sample material too, because the 
impedance mismatch affects the pressure of the reflected 
shock and therefore its transit time. It should be empha- 
sized that this 'ablator drive speed' does not represent 
any individual physical wave speed: it is a composite of 



FIG. 11: Pressure profiles at intervals of 1 ns for an Al sample 
coated with 20 fim of parylene-N ablator, flashed with 100 nm 
of Al, and driven with a laser pulse of lOPW/m^ for 3 ns. 
Before the start of the laser drive, the ablator extended from 
to 20 /im, and the sample from 20 to 45 /im. The profile at 1 ns 
shows the shock in the ablator. At 3 ns, the shock has been 
transmitted into the sample, almost doubling the pressure but 
exaggerating its deviation from fiatness. At 4 ns, the shock 
has almost reached the free surface of the sample, its peak has 
decayed and become less flat, and the ablation surface release 
is evident as a tilted step between 40 and 50 GPa on release. 
By 5 ns, the shock has reached the free surface of the sample 
and started to release. 
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FIG. 12: Pressure profiles at intervals of 1 ns for a Cu sample 
coated with 20 fim of parylene-N ablator, flashed with 100 nm 
of Al, and driven with a laser pulse of lOPW/m'^ for 3 ns. 
Before the start of the laser drive, the ablator extended from 
to 20 /im, and the sample from 20 to 45 /xm. The proflle 
at Ins shows the shock in the ablator. At 3ns, the shock 
has been transmitted into the sample, almost trebling the 
pressure but exaggerating its deviation from flatness. At 5 ns, 
the shock has almost reached the free surface of the sample, 
and the ablation surface release has caught up with the shock 
and made it triangular. 
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FIG. 13: Effect of ablator thickness on the pressure history 
applied to Al samples. Simulations for 10 /xm (solid) and 
20 fim (dashed) ablators, with a 3 ns drive pulse of 5, 10, and 
20PW/m^ 
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FIG. 14: Comparison between pressure induced by laser ab- 
lation in Al samples by a laser pulse 3 ns long of constant 
irradiance, with (solid) and without (dashed) a parylonc-N 
ablator. Irradianccs arc given in PW/m^. The pressure his- 
tory was calculated at the ablator-sample interface when the 
ablator was present, and at a Lagrangian position 1 /xm inside 
the sample when there was no ablator. 



the initial ablator shock speed, the speed of the reflected 
shock through the pre-compressed ablator, and the ab- 
lation rate. If the laser pulse is the optimum length for 
the ablator thickness, or longer, then the duration of the 
higher pressure drive experienced by the sample at the 
interface with the ablator can be related to the dura- 
tion of the laser pulse, as a load duration factor. If the 
laser pulse is longer than the optimum for a given abla- 
tor thickness, the principal pressure pulse applied to the 



FIG. 15: Comparison between pressure induced by laser ab- 
lation in Cu samples by a laser pulse 3 ns long of constant 
irradiance, with (solid) and without (dashed) a paxylene-N 
ablator. Irradiances are given in PW/m^. The pressure his- 
tory was calculated at the ablator-sample interface when the 
ablator was present, and at a Lagrangian position 1 /xm inside 
the sample when there was no ablator. 



sample does not become longer, but the sample pressure 
releases to the ablation pressure, and this lower pressure 
is sustained for longer. There is a minimum laser pulse 
duration for the peak pressure to reach the sample. If 
the laser pulse is longer than this, but shorter than the 
ablator thickness optimum, the duration of the pressure 
pulse at the sample varies proportionately from zero to 
the maximum. 

The initial mass density of the ablator is pQ. The 
ablation-driven shock travels at a speed Sa and com- 
presses the ablator to a mass density pa, ablating ma- 
terial at a rate ua with respect to the uncompressed ab- 
lator (the quantity calculated above) or Ua = UAPo/Pa 
of the compressed ablator. The shock reflected from the 
sample back through the ablator travels at a speed Sr 
and compresses the ablator further to a mass density pr- 
In this state, the speed of sound is c^: this is the speed 
of the head of the ablation surface rarefaction. For an 
ablator of initial, uncompressed thickness Iq, the transit 
time of the initial ablation shock is Ts — lo/sa'. this is the 
time that loading starts in the sample after the start of 
the laser pulse. The transit time of the reflected shock is 



Po/Pa - Ua/Sa 
) — • 



(9) 



The optimum laser pulse length is Ts+Tr, so the ablator 
drive speed is 



SaPa{Sr + Ua) 



(10) 



SaPO + SrPa 

The residual thickness of ablator when the reflected shock 
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reaches the ablation surface is 



Ir = lo [Po/Pa - Ua (Ts + Tr 



Pa 
Pr 



(11) 



SO the transit time of the ablation surface release is 



IpSr {poSg - PaUg) 



(12) 



The time for which the initial pressure is applied to the 
sample is thus 



loiPoSa - PaUa){PaSr + PrCr) 



[Sr + Ua) 

SO the load duration factor is 



Tl 



(poSa - PaUa){PaSr + PrCr) 
PrCriPQSa + PaSr) 



(13) 



(14) 



If the sound speed in the ablation shocked state is c^, 
traveling through the compressed ablator, the minimum 
pulse length for the peak pressure to reach the sample 
can be found by calculating the point in the ablator at 
which the head of a release wave at the end of the laser 
pulse would catch up with the ablation shock. The min- 
imum time is proportional to the sample thickness, so 
an effective 'catch-up speed' can be defined for a given 
ablation pressure: 



(15) 



PaCa - PoSa 



where Sa, Ca, and pa are all in the principal Hugoniot 
state for the ablator at the ablation pressure. For a given 
laser pulse duration r^, the maximum depth in the ab- 
lator before the ablation shock starts to decay is then 
UcTd- This calculation can also be used to calculate the 
maximum sample thickness for a supported shock in ex- 
periments where the sample itself is ablated ^2^. For an 
ablator of given thickness I, the minimum laser pulse du- 
ration Tj-^jj^ = l/uc- For laser pulse durations between 
this and Tmax = Ta + t^, the duration of the peak pres- 
sure pulse applied to the sample is 



fiTmax- 



Tmax — T„ 



(16) 



As a planar loading or unloading wave propagates 
through a sample, release waves from the edge will usu- 
ally erode the shocked region, [s^ The head of the release 
wave initiated by any lateral variation in the material or 
loading conditions is an expanding circle with respect to 
the shocked material (Fig. [TB)) : a cylinder or torus in 
two-dimensional plane or axisymmetric geometry respec- 
tively. The lateral variation in a projectile impact exper- 
iment is the edge of the projectile on the impact surface; 
for laser ablation it is the edge of the focal spot - though 
unlike a projectile impact, the pressure does not release 
to zero because the ablation of laterally-released mate- 
rial may induce a pressure almost as high as that from 



the shock-compressed ablator. The speed of sound in the 
shocked region is greater than the shock speed, so the 
release wave erodes the shock. The hydrodynamic flow 
is self-similar, so the erosion of the shock by the release 
wave can be characterized by the angle 4> at which the 
planar region of the shock is eroded. 4> is often assumed 
to be 45°, but this is not generally the case. As with 
the catch-up speed, these angles can be calculated given 
the sound speed on the principal Hugoniot Ca, by con- 
sidering the propagation of a disturbance from the edge 
of the planar shock region, traveling at Ca, across mate- 
rial moving with the particle speed Ua parallel with the 
shock, and the resulting speed of the disturbance as it 
moves across successive positions of the shock itself: 



tan( 



tan ( 



(Sq 



(17) 
(18) 



An equivalent analysis has been used previously to mea- 
sure the sound speed on the Hugoniot [l3|- For a circu- 
lar laser spot of diameter dj, the diameter of the region 
of the sample initially subjected to the maximum shock 
pressure is 



ds = di — 21 tan < 



(19) 



where / is the thickness of the ablator. Subsequent re- 
lease depends also on the sound and particle speed in 
the double-shocked ablator material. The same analysis, 
using the Hugoniot state in the sample, can be used to 
predict the lateral erosion of the shock in the sample it- 
self. The edge release angle is important when designing 
experiments to study phenomena at the leading edge of 
the shock. For experiments on release from the shocked 
state, including tensile damage, an important design pa- 
rameter is the scale rate at which lateral release propa- 
gates radially through the shocked material, /s„.[3l| 

The ablator drive speed, load duration factor, catch- 
up speed, edge angle, and lateral release ratio were calcu- 
lated using a generalized shock dynamics algorithm to de- 
termine the principal and secondary Hugoniot states [26j , 
using the same EOS as above (Figs [T71 to [^0]) . As was 
found previously for direct ablation of metals _25j, the 
catch-up speed exhibited a minimum value around the 
bulk modulus of the ablated material: at the minimum, 
the ablator thickness should be smallest for a given laser 
pulse duration. The ratio Ca/sa was greater than 3 for all 
relevant pressures, and rose past 10 for ablation pressures 
above 40 GPa: lateral release in the ablator may be a con- 
cern in experiments with thick ablators, though the pres- 
sure may not drop rapidly once release starts as it is sus- 
tained by the ablation. Radiation hydrodynamics sim- 
ulations were performed of the integrated system, with 
consistent results for the duration of the peak shock state 
in the sample. It should be noted that these quantities 
depend on material properties in states significantly off 
the principal shock Hugoniot of the ablator. The hydro- 
dynamic relations above are valid so long as radiation and 
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FIG. 16: Schematic of edge release eroding the shocked region. 
The lateral feature causing edge release is shown here as an 
interface (corner) with a material of lower shock impedance. 
The shock wave is shown at three positions: one just before 
interaction with the feature, where it is planar over the field of 
view, and two at different times after the interaction, showing 
the erosion of the shock by the rarefaction from the edge. 



charged particle transport are not significant within the 
compressed ablator or the sample. The radiation hydro- 
dynamics simulations indicated that radiation transport 
is unlikely to be significant in this regime; they did not 
include electron transport. For ablatively driven experi- 
ments on a range of materials at a range of shock states, a 
wide range of off-Hugoniot states may be involved. These 
properties are often not known at all accurately from ex- 
perimental measurements. It would be more efficient to 
use experimental measurements to validate and adjust 
wide-ranging theoretical models, rather than attempting 
to measure the properties of the ablator so generally. Ex- 
periments should include the principal Hugoniot of the 
ablator, and also off-Hugoniot states induced by, for ex- 
ample, reflected shocks from high-density samples. These 
data could be acquired using laser loading experiments 
with velocimctry measurements. Transparent samples 
would allow the velocity history to be measured at the 
interface. Whether transparent or opaque, the samples 
should be thick enough to prevent release waves from 
their free surface from perturbing the states inside the 
ablator, which would complicate the experiments. Mea- 
surements would include the amplitude and duration of 
the shock in the sample. 



VI. CONCLUSIONS 

The sensitivity of ablation pressure in CH-based plastic 
ablators to uncertainties in equation of state or composi- 
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FIG. 17: Example calculations of 'ablator drive speed' relat- 
ing the optimum ablator thickness and laser pulse duration, 
and load duration factor relating the laser pulse duration and 
pressure pulse duration on the sample, for a parylene-N abla- 
tor and Al or Cu samples. 
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FIG. 18: Example calculations of 'catch-up speed' relating 
the maximum ablator thickness and laser pulse duration, for 
a parylene-N ablator. 



tion was predicted to be 5-10%, which is small compared 
with typical uncertainties in laser irradiance. A laser 
pulse of constant irradiance was predicted to induce an 
ablation pressure that decreased by several percent per 
nanosecond, but was significantly more constant than the 
pressure history induced by direct ablation of metal sam- 
ples. A simple relationship between laser irradiance and 
shock pressure was deduced from the simulations, valid 
for irradiances from 0. 1-100 PW/m^. Al flashing, used to 
prevent early-time shine-through, was predicted to make 
at most a few percent difference in loading history once 
dielectric breakdown had occurred, increasing the pres- 
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FIG. 19: Shock erosion angle, describing the angle at which 
the lateral release wave propagates across the shock, for a 
parylene-N ablator. 
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FIG. 20: Ratio between bulk sound speed and shock speed, 
describing the rate at which the lateral release wave propa- 
gates radially through the shocked materal, for a parylene-N 
ablator. 



sure and flattening the pressure history at lower irradi- 

anccs. It was demonstrated that the irradiancc history 
can be adjusted to produce a more constant pressure his- 
tory 

The impedance mismatch between plastic ablators and 
metal samples induces a stronger transmitted shock in 
the sample, by a factor of around 1.8 for Al and 2.5 for 
Cu. The shock reflected into the ablator interacts with 
the ablation surface to produce a stepped pressure his- 
tory in the sample, the pressure reducing to closer to the 
ablation pressure. This composite loading history should 
be taken into account when designing and interpreting 
material dynamics experiments. The shock pressure is 
applied to the sample for a significantly shorter time than 
the laser pulse duration, because the successive waves 
in the ablator are faster than the initial ablation-driven 
shock. 

A compact method was found for representing the op- 
timum ablator thickness for a given laser pulse dura- 
tion (or vice versa), as an effective speed. A similar 
method was found to represent the relationship between 
the laser pulse duration and the duration of the initial, 
peak pressure pulse applied to the sample. These rela- 
tions are helpful in the design and interpretation of ab- 
latively driven shock experiments. 
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